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Abstract 

In the framework of perturbation theory the reality of the perturbed eigenvalues 
of a class of "PTsymmetric Hamiltonians is proved using stability techniques. 
We apply this method to "PTsymmetric unperturbed Hamiltonians perturbed by 
"PTsymmetric additional interactions. 
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§ ■ 1 Introduction 

o 

Perturbation theory has played a very important role in the past ([1]) in the study 
of non Hermitian Hamiltonians with VT symmetry ([2]). In [3, 4, 5] it was applied 
starting from a selfadjoint Hamiltonian in order to investigate the perturbation of its 
spectrum by a VT symmetric interaction. In the present paper we try to present 
more general results concerning in general the perturbation of a non selfadjoint VT 
^ ■ symmetric Hamiltonian: 

c3 ; H ^ H*, HVT = VTH . 

This extension of perturbation theory will be shown to have non trivial aspects even 
when one restricts oneself to unperturbed discrete simple levels with real energies. 
Our aim is to provide a consistent and selfcontained framework where well established 
results and new developments can be discussed coherently. The main original aspects 
that we present are: 

1. extension of perturbation theory for VT symmetric Hamiltonians in order to 
prove the reality of the spectrum of a class of VT symmetric operators perturbed 
by a VT symmetric interaction; 

2. the interaction does not need to be bounded relative to the unperturbed Hamil- 
tonian; 
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3. the results concerning the reality of the eigenvalues are achieved by using the 
stability theory developed by Hunziker and Vock (HV-stability theory) in [6]. 

The paper is organized as follows. In Section 2 we give a general presentation of 
perturbation theory for non selfadjoint VT symmetric operators and a review of the 
stability theory for eigenvalues. In Section 3 we give the technical results leaving the 
proof to Section 4. Some open problems and further perspectives are outlined in the 
Conclusions (Section 5). 

2 General formalism in perturbation theory and re- 
view of results 

An operator is called VT symmetric if it is invariant under the combined action of a 
reflection operator V and the antilinear complex conjugation operator T . One basic 
issue is to prove the reality of the spectrum. Results in this direction have been obtained 
in [7, 8] by means of ODE techniques, and they have been recently extended in [3, 4, 5] 
(see also [9] for a brief review) in the framework of perturbation theory by perturbing 
selfadjoint Hamiltonians. The main goal of this paper is to extend part of these results 
to the case of VT symmetric Hamiltonians obtained by perturbing a VT symmetric 
operator, not necessarily selfadjoint. Thus, for the convenience of the reader, we first 
recall the results of [3, 4, 5], concerning VT symmetric operators in a Hilbert space 
(H, < u, v >) of the form 

H g = H + igW, 9 6R, 

where Hq is selfadjoint, with domain D C 7i, and W is a symmetric operator, relatively 
bounded with respect to H . Moreover there exist a unitary involution V : 7i — > 7i 
and an antilinear involution T : 7i — > H, both mapping D to D, such that 

H V = VH , VW = -WV, TH = H T, TW = WT . (2.1) 

Then H g is VT symmetric, i.e. VTH g = H g VT, for g e R. An example is provided 
by VT symmetric Schrodinger operators in L 2 (R d ), d > 1, where H is the selfadjoint 
realization of —A + V, T is complex conjugation, i.e. (Tip)(x) = i/)(x), and V is the 
parity operation defined by 

(Vi>){x) = (-l) jd x d ) , V e L 2 (R d ) (2.2) 

where jk = 0, 1, and jk = 1 for at least one 1 < k < d. Here —A denotes the d- 
dimensional Laplace operator; V and W are real valued functions, P-even and V-odd 
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respectively: VV = V,VW = —W. The following theorem, proved in [3, 4], provides 
a result in the case of bounded perturbation. 

Theorem 2.1 Let H be a self adjoint operator in H and W a symmetric operator in H, 
satisfying the above assumption (2.1). Assume further that H is bounded below, W is 
bounded and that the spectrum of H is discrete. Let a (Ho) = {Ej : j — 0, 1, . . .} denote 
the incresing sequence of distinct eigenvalues of Ho- Finally, let 5 := infj> [Ej + i — 
Ej}/2 and assume that 5 > 0. Then the following results hold: 

(i) if for each degenerate eigenvalue of Hq the corresponding eigenvectors have the 
same V -parity, i.e. they are either all V-even or all V-odd, then cr(H g ) C R, if 
geR, \g\<5/\\W\\; 

(ii) if H has an eigenvalue E with multiplicity 2 whose corresponding eigenvectors 
have opposite V -parity, i.e. one is V-even and the other one is V-odd, then H g 
has a pair of non real complex conjugate eigenvalues near E for \g\ small, g G R. 

The proof of the reality of the spectrum under suitable conditions has been extended 
to the case of relatively bounded perturbation in [5] and is recalled in the following 

Theorem 2.2 Let 

d 2 

H g = -— + V(x) + igW(x) in L 2 (R) , 

where V(x) is a real-valued even polynomial of degree 21, with lnm-ri^oo V (x) = +oo 
and W(x) is a real-valued odd polynomial of degree 2r — 1, with I > 2r. Then there 
exists go > such that o~(H g ) C R for \g\ < g . 

Remark 2.3 The above theorems enlarge the class of 'PTsymmetric Hamiltonians 
with real spectrum provided in [7, 8]. 

The question that we want to address in this paper is the folllowing. Let now H 
be VT symmetric (not necessarily selfadjoint) with real spectrum and let W be VT 
symmetric as well. We will provide conditions on H and W in order to guarantee 
that the perturbed eigenvalues of H(e) := H + eW, e real, are real for |e| small. 
Therefore we shall not prove that the whole spectrum of H(e) is real, but that at least 
the eigenvalues generated by the unperturbed real ones stay real. First of all notice 
that if H has degenerate eigenvalues there may be problems even in the selfadjoint 
stated in Theorem 2.1(h). In general we have 

m g (X) < m a (X) (2.3) 
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where 

m g (\) = dim{w : (H - X)u = 0} (2.4) 
is the geometric multiplicity of A G C and 

m a (X) = dim{w : (H - X) n u = , for some n G N} (2.5) 

is the algebraic multiplicity of A. Notice that in the finite dimensional case (i.e. H and 
W are matrices) m a (X) is the multiplicity of A as a root of the characteristic polynomial 
of Hq. If Hq is selfadjoint then m g (X) = m a (X), for all A G C. However in the general 
(non selfadjoint) PTsymmetric case it is not enough to assume that m g (X) = 1 in 
order to guarantee that the eigenvalues are simple, i.e. that there is no degeneracy, i.e. 
that m a (X) = 1. In fact there may be "exceptional points", caused by the presence of 
Jordan blocks. 

Example 2.4 Set 

(26) 

and let T be complex conjugation. Then Hq is VT symmetric. Its characteristic 
polynomial det(if — XI) = A 2 has just one root A = with m 9 (A ) = 1 and m a (A ) = 
2. If we take 

"-(?'.) 

as a perturbation, then H(e) = H + eW has a pair of non-real complex conjugate 



eigenvalues A(e) = ±iy e(e + 2) for e > 0. 

To avoid such difficulties we analize here only the non-degenerate case, i.e. we assume 
that 

m a {X) = 1 (2.7) 

for all eigenvalues A of H . The Hamiltonians H that we perturb are one-dimensional 
PTsymmetric Schrodinger operators with real simple spectrum, for example the Hamil- 
tonians H g provided by the above Theorems 2.1 and 2.2 for \g\ small. 

Remark 2.5 Perturbation theory for matrices in the non-degenerate case provides a 
straightforward result because: 

(1) the stability of the unperturbed eigenvalues (see Definition 2.6 below) is guar- 
anteed by the boundedness of the perturbation W (any matrix W is a bounded 
operator). This implies that near any unperturbed eigenvalue A of Hq there is 
one and only one eigenvalue A(e) of H(e) for |e| small and lim e _^ A(e) = A. Thus 
m a (X(e)) = 1; 
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(2) the eigenvalues of PTsymmetric operators come in pairs of complex conjugates, 
i.e. if A(e) is an eigenvalue of H(e), then A(e) is an eigenvalue of H(e) too. 

Therefore, combining (1) and (2), A(e) G R. 

We will provide new results on the reality of the perturbed eigenvalues (A(e)) in the 
case of unbounded (not necessarily relatively bounded) perturbation. As observed in 
Remark 2.5, a crucial issue is the stability of the unperturbed eigenvalues, i.e. the 
problem is reduced to a stability result for the eigenvalues of Hq. Such result is imme- 
diate if the perturbation W is bounded relative to Hq (see e.g. [10, 11]). Indeed, in 
this case the Rayleigh-Schrodinger perturbation expansion (RSPE) is convergent, for 
|e| small, to the (real) perturbed eigenvalues. However there can be stability even if 
the RSPE does not exist. We provide here a brief reminder on stability (see [10] and 
[6])- 

Definition 2.6 A discrete eigenvalue A of H is stable with respect to (w.r.t.) the 
family H(e) = H + eW if 

(i) for any small r > 0, 

r r = {z :| z - A |= r} C p(H(e)) , as e -> 
where p{H{e)) := C — a(H(e)) is the resolvent set of H(e); 

(ii) lim £ _ || P(e) - P(0) ||= where 

P(e) = (2m)- 1 £(z- H(e))- 1 dz 

is the spectral projection of H(e) corresponding to the part of the spectrum en- 
closed in T r , and H(0) := H Q . 

Since (ii) implies that 

dimP(e) = dimP(O) (= m (A)) (2.8) 

for |e| small, A is the limit of a group of perturbed eigenvalues with the same total 
algebraic multiplicity (see [10]). If W is bounded relative to H there exist a, b > 
such that 

|| Wu \\< b || H u II +a || u || , VueD(H ). (2.9) 

Then 

|| (zo-Hie))' 1 -(zo-Ho)' 1 ||->0, as e -> (2.10) 
for some z a(H ) and this implies (ii). 
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Remark 2.7 Eigenvalues may be stable even if (2.10) fails. As an example let H(e) = 
p 2 +x 2 +ex A , e > 0, in 7i = L 2 (TV) be the Hamiltonian of the even anharmonic oscillator. 
Here p 2 = —d 2 /dx 2 . Then the eigenvalues of the harmonic oscillator H(0) are stable 
w.r.t. H(e), e > 0, in spite of the fact that (2.10) fails and the RSPE is divergent. In 
other words, the continuity of the eigenvalues at e = holds although analitycity fails. 
In this case we only have strong resolvent convergence, i.e.: 

\im(z -H(e))- 1 u=(z -H(0))- 1 u, Vw G H (2.11) 

for some zq a(H(0)). Then (2.11) yields the strong convergence of the projections 
P(e), i.e. 

lim P(e)u = P(0)u , Vw G H (2.12) 

which implies 

dimP(e) > dimP(O) , as e -> . (2.13) 

So H(e) may have more eigenvalues than H(0) in the circle T r . This happens for 
instance for the double-well operator 

H(e) =p 2 + x 2 (l-ex) 2 , in L 2 (R). (2.14) 

In fact near any eigenvalue of the harmonic oscillator P^(0) there are two eigenvalues 
of H(e) for |e| small. 

3 Statement of the results 

Let H(e) = p 2 + V + eW, e > 0, denote the closed operator in L 2 (K) with C£°(R) as a 
core, defined by 

H{e)u = — u" + Vu + eWu , Vw G D(H(e)) , (3.1) 

where V — V + + iV-, W = W + + iW- and V+, V-, W + , W- are real valued functions 
in L^ C (R). Moreover V + , W + are P-even and bounded below, V_, W- are P-odd and 

lim \V(x) + eW{x)\ = +oo. (3.2) 

| x | — *oc 

Here (Vu)(x) = u(—x),\/u G L 2 (R). Let us further assume that the spectrum of H(e), 
denoted a(H(e)), is discrete for e G [0,e ], i.e. it consists of a sequence of isolated 
eigenvalues with finite algebraic multiplicity: 

a(H(e)) = {E J (e):j = 0,l,...}. (3.3) 
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For the unperturbed eigenvalues (Ej(0)) we will adopt the simplified notation Ej : = 
Ej^J = 0,1,... 



Theorem 3.1 Under the above assumptions the following statements hold: 

(1) each eigenvalue Ej of H(0) is stable w.r.t. the family H(e),e > 0. In particular, 
if Ej is simple, i.e. m a (Ej) = 1, and real there exists ej > such that for 
< e < €j, H(e) has exactly one eigenvalue Ej(e) close to Ej.- 

lim Ea (e) = Ej 

and Ej(e) is real; 

(2) there are no "dying eigenvalues" , i.e. if E(e) G a(H(e)) and lim e ^o E{e) = E, 
then E is an eigenvalue of H(0). 

In the following examples all the above conditions are satisfied. Moreover, the eigen- 
values of H(0) are real and simple. 

Example 3.2 The unperturbed Hamiltonian H(0) can be any of the operators H g 
satisfying the above Theorems 2.1 and 2.2, for instance H(0) = p 2 + x 2n + igsinx, 
H(0) = p 2 + x 2n + igx/ (x 2 + 1), H(0) = p 2 + x 21 + igx 2q ~ x , with \g\ small, n, I, q G N 
and / > 2q. Another example is provided by H(0) = p 2 + ix 2k+l , k G N. In all 
these cases the perturbation W can be taken of the form W = W + + iW- where 
W + = exp (x 2 ), or W + is an even polynomial function diverging positively at infinity, 
and W- is an odd polynomial function bounded relative to W + . If -f^(O) = p 2 + ix 2k+1 , 
k G N, the polynomial W- must diverge positively at +oo in order to guarantee that 
condition (3.2) is satisfied. 

Corollary 3.3 Assume that the eigenvalues Ej,j G N 7 of H(0) are all simple and 
real. Then all the perturbed eigenvalues Ej(e),0 < e < €j, of H(e) are real and simple. 
Moreover (non-real) complex eigenvalues of H(e) cannot accumulate at finite points but 
only at infinity. 

Remark 3.4 1. The question of the reality of the whole spectrum of H(e) is still 
open, because there may be (possibly complex) eigenvalues of H(e) diverging to 
infinity as e — > 0. 

2. Hamiltonians satisfying Corollary 3.3 are provided in the above Example 3.2. 
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If the perturbation W is bounded relative to H(0) then the RSPE near any Ej con- 
verges. We have therefore the following 

Corollary 3.5 If the unperturbed eigenvalues Ej are simple and real and W is bounded 

relative to H(0), then the RSPE near Ej is real for all j — 0, 1, More precisely, for 

< e < ej 

oo 

£i(f) = Ea/, a = E J (3.4) 

?1=0 

and a n G R, Vn G N. 

Remark 3.6 If W is not relatively bounded w.r.t. H(0), the RSPE although divergent 
can be Borel summable to Ej(e). In this case the coefficients of the RSPE are real. This 
happens for instance for the Hamiltonians H (0) and the perturbations W of polynomial 
type described in Example 3.2. 

Remark 3.7 The stability theory developed in [6] allows one to prove a result similar 
to that stated in Theorem 3.1 also in presence of essential spectrum provided that 

dist(£y, a ess (H(e))) > c> , as e -> 0. (3.5) 

where a ess (H(e)) is the complement in a(H(e)) of the discrete spectrum of H(e). 

4 Proof of Theorem 3.1 

Although it is a straightforward application of the HV-stability theory developed in [6], 
for a pedagogical purpose and for the convenience of the reader we describe here the 
main steps. The only " continuity condition" required by the HV-theory and guaranteed 
by the assumptions of the theorem is: 

KmH(e)u = H(0)u, G C °°(R). (4.1) 

Condition (4.1) implies the strong convergence of the resolvents (2.11) which yields 
(2.13) and is not enough, as already remarked, to ensure stability. One key ingredient 
is the numerical range of H(e): 

N(e) := {< u, H{e)u >: u G D(H(e)), \\u\\ = 1} (4.2) 

which contains the eigenvalues of H(e) and in the present case is contained in the right 
half-plane: 

a(H(e)) C N(e) d TZ+ := {z : Re z > 0}. (4.3) 
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Next we need to introduce the set of uniform boundedness of the resolvents: 

V := {z : (z — H(e))~ l exists and is uniformly bounded as e — > 0}. 

The sets V and N(e) are closely related to each other (see e.g. [10]), since C— N(e) C V. 
In particular for z G 1Z_ := {z : Rez < 0}, we have 

\\(z - Hie))- 1 ]] < dist(z,7V(e)) _1 < \Rez\- 1 , Ve > 0. (4.4) 

Thus z G T>. The HV-theory allows us to prove that T> is much wider than 1Z-. Indeed 
we can prove that the following alternative holds: 

(1') if E G <t(H{0)) then E is stable w.r.t. H{e), e > 0; 
(2') if E i cr(H(0)) then E G V. 

Thus, T> coincides with the complement of the spectrum of H(0). Notice that (2') 
implies statement (2) of the theorem. Moreover a first step in the proof of (1') consists 
in showing that the circle r r is contained in £>, and this guarantees the possibility to 
construct the projection -P(e) for all e sufficiently small (see the above Definition 2.6). 
Finally, a third key ingredient, on which the proofs of both (1') and (2') are based, 
is represented by the so called "characteristic sequences" of (z — H(e)) or "Weyl-type 
sequences", i.e. sequences (e n ,u n ) such that: 

e n > 0, u n G D(H(e n )), \\u n \\ = 1, u n A 0, \\(z - H{e n ))u n \\ -> 0. (4.5) 

Now it is easy to check that for a suitable constant a > we have 

< u, p 2 u >< a(Re < u, H(e)u > + < u, u >) , (4.6) 

Vw G D(H(e)), Ve G [0,e ]. By making use of (4.6) one can prove that a sequence 
of type (4.5) generates another characteristic sequence (e n ,v n ) which is "supported at 
infinity", i.e. such that 

v n (x) = , for \x\ < n , Vn G N. (4.7) 

Since limui^oo |y(x)| = +oo by the above assumption (3.2), we first assume that 
limia-i^oo V + (x) = +oo. Now the idea is to prove (2'), i.e. that a(H(0)) U T> = C, 
by contradiction. More precisely, one can show that if z ^ a(H(0)) U V, then a 
characteristic sequence (e n ,v n ) of (z — H(e)) exists and it satisfies both (4.5) and (4.7). 
But this contradicts the fact that for all z G C 

lim d n (z, e) = +oo (4.8) 

e->0 
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where 



d n (z,e) = w£{\\(z - H(e)) : u G D(H(e)) , \\u\\ = l,u(x) = Ofor \x\ < n} . (4.9) 
Indeed we have 

dn(z,e) > dist(z,N n (e)) (4.10) 

where 

N n (e) = {< u, H(e)u >: u G D(H(e)), \\u\\ = 1 , u(x) = for \x\ < n} (4.11) 

is the so called "numerical range at infinity". Now, for all u G D(H(e)) such that 
= 1 and u(x) = for \x\ < n, we have 

\z— < u, H(e)u > | > | < u, H(e)u > \ — \z\ 
> | Re < u, H(e)u > | - \z\ > < u, V + u > —\z\ — c (4-12) 

as e — > 0, for some constant c > 0. Thus, (4.8) follows from (4.10) and the fact that 
lim^i^oo V+(x) = +oo by assumption. Hence (2') is proved. Now let E G <r(if(0)). In 
order to prove that E is stable w.r.t. H(e), e > 0, we recall that it is enough to prove 
(2.8), i.e. in view of (2.13), that 

dimP(e) < dimP(0) , as e -> . (4.13) 

Again the proof is by contradiction. In fact if we assume that dimP(e„) > dimP(0), 
for e n — > 0, a characteristic sequence of (E — H(e)) which satisfies (4.7) can be found, 
and exactly as before this contradicts (4.8). Thus the theorem is proved under the 
assumption that lim\ x \^, 00 V + (x) = +oo. If this assumption is not satisfied then it 
follows from (3.2) that lim^i^ |Vl(x) + eW^_(:r)| = +oo. Without loss of generality, 
since V- and W- are P-odd, we may assume that \im x ^ +00 (V-(x) +eW-(x)) = +oo and 

e— >0 

lim a: ^-o (Vl(a;) + eW-(x)) = — oo. This immediately implies that there exists c > 
such that 

lim d n (z, e) = , Vz G 1Z+ , Re z > c 

and therefore (4.8) cannot be used to generate a contradiction with (4.5) and (4.7). 
Nevertheless the problem is overcome as follows (see also [12] where an analogous prob- 
lem was handled in a similar fashion). First of all one can prove that any characteristic 
sequence (e n ,u n ) generates another characteristic sequence (e n ,v n ) which is supported 
either at +oo or at — oo, i.e. such that either 

v n (x) = , for x < n , Vn G N (4.14) 
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or 

v n (x) = 0, for x > -n, Vn G N. (4.15) 

Next we introduce the numerical range at +00, iV+(e), and the numerical range at 
—00, N~(e), defined by (4.11) where the condition u(x) = for \x\ < n is replaced by 
u(x) = for x < n, and u(x) — for x > — n respectively. More precisely: 

N+(e) = {< u, H(e)u >: u G D(H(e)), \\u\\ = 1 , u(x) = for x < n} 

and 

^C( e ) = {< u,H(e)u >:ue D(H(e)), \\u\\ = l,u(x) = Ofor x > -n} . 
Then for all z G C we have 

lim d±(z,e) = +00 (4.16) 

where 

d+(z, e) = inf{||(*- if (e)) : u G D(H(e)) ,\\u\\ = l,u(x) =0fovx<n}. (4.17) 

and similarly for d~(z,e). Then the existence of a characteristic sequence satisfying 
either (4.14) or (4.15) contradicts (4.16) and this concludes the proof of the theorem. 

5 Conclusions 

The reality of the spectrum of a PTsymmetric Hamiltonian has been proved for a 
class of polynomial Hamiltonians and for solvable classes of potentials generated by 
solvable selfadjoint problems by a complex coordinate shift. By a suitable use of the 
perturbation theory we have summarized in this paper that it is conceivable that the 
class of PTsymmetric Hamiltonians with real spectrum can be considerably enlarged 
in so far as one starts from a PTsymmetric Hamiltonian with real spectrum and adds 
a VT symmetric perturbation. This analysis does not pretend to be exhaustive in the 
sense that the perturbed eigenvalues of H (e) are real but there is no guarantee that 
the whole spectrum is real. Therefore the connection with pseudo-Hermiticity deserves 
further investigation. Moreover the pseudo-Hermiticity condition H = i]H*i]^ 1 (see e.g. 
[13, 14]) of the unperturbed Hamiltonian could generate a perturbative approach to 
pseudo-Hermiticity for larger classes of Hamiltonians and of the operator 77. Another 
open question is a detailed analysis of the degenerate case, i.e. the case when the 
unperturbed eigenvalues are diabolic points ([15]) or exceptional points ([16, 17]); the 
latter are typical of non diagonalizable Hamiltonians. 



11 



Acknowledgements 

We wish to thank D. Heiss, H. Geyer and M. Znojil for the warm hospitality at the 
workshop and for providing the opportunity to publish the contributions. We also 
thank P. Dorey for useful discussions and suggestions. 

References 

[I] Caliceti E, Graffi S and Maioli M 1980 Commun. Math. Phys. 75 51 

[2] Bender C M and Boettcher S 1998 Phys. Rev. Lett. 80 5243; Cannata F, Junker G 
and Trost J 1998 Phys. Lett. A 246 219; Andrianov A A, Cannata F, Dedonder 
J P and Ioffe M V 1999 Int. J. Mod. Phys. A 14 2675; Bender C M, Boettcher 
S and Meisinger P 1999 J. Math. Phys. 40 2201; Bender C M, Milton K.A. and 
Savage V M 2000 Phys. Rev. D 62 085001; Bender C M, Brody D C and Jones 
H F 2004 Phys. Rev. Lett. 93 251601; Znojil M 1999 Phys. Lett. A 259 220 

[3] Caliceti E, Graffi S and Sj6strand J 2005 J. Phys. A: Math. Gen. 38 185 

[4] Caliceti E, Graffi S and Sjostrand J: in preparation 

[5] Caliceti E and Graffi S 2005 J. Nonlinear Math. Phys. 12, Supplement 1, 138 

[6] Hunziker W and Vock E 1982 Commun. Mathy. Phys. 83 281 

[7] Dorey P, Dunning C and Tateo R 2001 J. Phys. A: Matyh. Gen. 34 L391 

[8] Shin K C 2002 Commun. Math. Phys. 229 543 

[9] Caliceti E 2004 Czech. J. Phys. 54, n 10, 1065 

[10] Kato T 1976 Perturbation Theory for Linear Operators, 2nd Edition (Berlin: 
Springer) 

[II] Reed M and Simon B 1978 Methods of Modern Mathematical Physics vol IV 
(New York: Academic) 

[12] Caliceti E and Maioli M 1983 Ann. Inst. H. Poincare A 38 175 



12 



Mostafazadeh A 2002 J. Math. Phys. 43 205; Mostafazadeh A 2003 J. Phys. A: 
Math. Gen. 36 7081; Jones H. and Mateo J. 2005 Czech J. Phys. 55, n. 9, 1117; 
Mostafazadeh A 2005 Czech J. Phys. 55, n. 9, 1157 

Figueira de Morisson Faria C and Fring A Preprint quant-ph 0604014; Bagchi B, 
Banerjee A and Quesne C Preprint quant-ph 0606012 

Kirillov O and Giinther U 2006 Preprint math-ph 0602013 

Seyranian P A, Kirillov O N and Mailybaev A A 2005 J. Phys. A: Math. Gen. 
38 1723 

Weigert S 2005 Czech J. Phys. 55, n. 9, 1183 



13 



